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Quantum phase transition from the Neel to the dimer states in an antiferromagnetic(AF) Heisenberg 
model on square lattice is studied. We introduce a control parameter a for the exchange coupling 
which connects the Neel (a = 0) and the dimer (a = 1) states. We employ the CP representation 
of the s = i spin operator and integrate out the half of the CP 1 variables at odd sites to obtain a 
CP 1 nonlinear a model. The effective coupling constant is a function of a and at a = the CP 1 
model is in the ordered phase which corresponds to the Neel state of the AF Heisenberg model. A 
phase transition to the dimer state occurs at a certain critical value of ac as a increases. In the 
Neel state, the dynamical composite U(l) gauge field in the CP 1 model is in a Higgs phase and 
low-energy excitations are gapless spin wave. In the dimer phase, a confinement phase of the gauge 
theory with s = 1 excitations is realized. For the critical point, we argue that a deconfinement 
phase, which is similar to the Coulomb phase in 3 spatial dimensions, is realized and s = \ spinons 
appear as low-energy excitations. 



Quantum phase transition(QPT) is one of the most 
interesting problem in these days. 1 It is often argued 
that the simple Ginzburg-Landau theory does not ap- 
ply to certain class of the QPT's. In this paper we shall 
study s = \ antiferromagnetic(AF) Heisenberg model on 
2-dimensional square lattice with nonuniform exchange 
couplings, 
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where x denotes site of the spatial lattice, j is the direc- 
tion index (J = 1,2) and S x is the spin operator at site 
x. We rename the even lattice sites x = (o, i) where o 
denotes odd site and the index i = 1, 2, 3 and 4 specifies 
its four nearest-neirfiborfNN") even sites (see Fie.l). 



where < a < 1 is a control parameter which connects 
the uniform Heisenberg model to the dimer model. 

It is not so difficult to derive the CP 1 field-theory 
model 2 ' 3 from Eq.(l) 4 . The spin operator S x can be ex- 
pressed in terms of the CP 1 variable z x = {z\,z^) t as 
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where a are the Pauli matrices and the CP 1 constraint 
S a =i 2 \ z x\ 2 = 1 guarantees the magnitude of the local- 
ized spin as \. 

From our assignment of J Q i (2), it is obvious that J Q \ 
is larger than the others. We use the path-integral for- 
malism and parameterize the CP 1 variable z Q by refering 
to Z i = z e , 
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FIG. 1. Numbering and dimer picture:cross symbol is odd 
site, and spot symbol is even site. Solid line means that its 
correlation is stronger than dotted lines. 



The exchange couplings J X j = J i are position depen- 
dent and we explicitly consider the following case which 
corresponds to the dimer configuration, 
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where p is a parameter, e is a phase factor and 
z e = io-2Z*. At vanishing temperature (T), spins tend to 
point antiparallel their NN spins, and then the parame- 
ter p can be treated as a small parameter. We expand 
y/l — \p \ 2 ~ 1 — \ \p \ 2 + ■ ■ ■ and retain only terms up to 
quadratic of p . Then we perform the Gaussian integra- 
tion of p 's to obtain an effective model of z e 's for which 
smooth configurations dominate at T — 0. 

Calculation is rather long but straightforward 4 and 
we obtain an effective field theory of the AF Heisenberg 
model under study, 

£cp= [l^l'+^kl 2 --^)] +C B , (5) 
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where D^Zx = (<9 M + iA^)z x , = iz^d^z is the U(l) 
gauge held, a x is the Lagrange multiplier for the CP 1 
constraint and Cb is the Berry phase term. We have 
rescaled the imaginary time r and x as taking the con- 
tinuum limit (the spatial coordinate x in Eq.(5) denotes 
the even site of the original lattice). The effective cou- 
pling constant / e g and the "speed of light" c, which is 
often set unity, are explicitly given as 
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where a is the lattice spacing of the original lattice and 

/AFH - 2V2a and c AF H = ^r 1 - 

Here we should comment on the Berry phase terms, 
Cb, which appear in the CP 1 nonlincar-er model repre- 
sentation of the AF Heisenberg model. First we consider 
AF spin chains in one dimension in order to obtain impor- 
tant insight for the effect of the Berry phase terms. For 
the uniform AF chains, the Berry phase reduces to the 
#-term (i6 J drdxF TX ) with 9 — it. The gauge dynam- 
ics crucially depends on the value of 9 and only the case 
of 9 = 7r (mod 2ir), the deconhnement phase is realized. 
Topologically nontrivial configurations are suppressed by 
the #-term when 9 = n. 

On the other hand for the bond-alternating(BA) AF 
Heisenberg chain with J x = J + (—1) X AJ, the value of 
the 9 parameter is 9 = jq^j Tr, and therefore the above 
suppression by the #-term does not occur and gauge dy- 
namics is essentially the same with the 9 = case. Then 
the confinement of the spinon z x occurs, and low-energy 
excitations are gapful s = 1 excitations n x = ^z x az x 
even for a very small value of AJ. This is nothing but 
the spin-Peierls transition. 

For the 2-dimcnsional uniform AF Heisenberg model, 
the Berry phase term similarly appears as Cb = 
|Ei £ i^i where e x = 1(— 1) for even (odd) site x (of 
the lattice with the lattice spacing y/2a for the CP 1 
model (5)) and A x is the area enclosed by the curve given 
by the time evolution of n x (T). It can be proved that it 
gives only vanishing contribution for an arbitrary smooth 
configuration. However for singular configurations like 
instantons (monopolcs), the Berry phase is nonvanishing 
and it gives terms like Cb = *f Eti EiCKj where 
£* are 0, 1, 2,3 (i = 1, 2, 3, 4), respectively and m x is the 
instanton number. 5 ' 6 (The index i refers to the four dual 
sublattices of the original lattice.) From the above Cb, it 
is obvious that the cancellation mechanism of the instan- 
ton contribution occurs by the destructive interference 
unless m x = 4 (mod 4) . 

For the nonuniform cases which are considered in this 
paper, the coefficient of the Berry phase term becomes ir- 
rational as in the case of the BA spin chains and therefore 
cancellation of instanton effect does not occur or at least 
weakens. Then we can expect that the (global) phase 



structure of the effective gauge model for the nonuni- 
form AF Heisenberg model is the same with that of the 
CP 1 model without the Berry phase terms. More com- 
ments will be given after showing results of nonpertur- 
bative studies on the gauge dynamics. 

It is straightforward to obtain the effective potential of 
a = ( a x) and z 2 = (z x ) by integrating out z x , 

V cS = o-(\z 2 \ 2 -^)+^-[(a+A 2 )^+A^-A 3 -«V^ 

where the cutoff A = The effective potential V c r 

indicates that there exists a critical coupling fc = x- 
The existence of the phase transition has been verified by 
the numerical calculation of the equivalent 0(3) nonlin- 
ear a model in (2 + 1) dimensions 7 . In the weak-coupling 
region / c gf < fc, the spontaneous symmetry breaking oc- 
curs and (z x ) 7^ 0. As a result, the Higgs phase is realized 
in the gauge-theory terminology. Low-energy excitations 
are gapless spin waves which are described by z x . In 
the strong-coupling phase / e g- > fc, (o~ x ) ^ whereas 
(z x ) = 0. Local Maxwell terms appear in the effective 
action of the gauge field A^, and the confinement phase 
is realized. Low-energy excitations are s = 1 compos- 
ite of the spinons which correspond to n x with (mass) 2 

OC (<7 X ). 

The effective coupling / c g in (6) first decreases as a in- 
creases but above certain value of a it starts to increase 
and goes to infinity at the dimer limit a = 1. In the 
uniform case a = 0, f e s(a = 0) < fc and this means 
that the ordered Neel state is realized at the vanishing T 
in the AF Heisenberg model in two spatial dimensions as 
it is now widely believed. The behavior of the effective 
coupling / e ff(a) shows the existence of a critical value 
ac at which the phase transition occurs. This result in- 
dicates that the strong-coupling phase of the CP 1 model 
(5) corresponds to the dimer phase in which the ground 
state is nothing but spin-singlet pairs formed by the al- 
ternative strong exchange couplings and excitations have 
s = 1. In fact this Neel-dimer transition was observed by 
the numerical calculations some years ago 8 . 

Hereafter we are interested in the critical point at 
feS — fc which separates the Neel and dimer phases. 
In order to investigate that "phase" , study on the gauge 
dynamics is required. At fc, (<J x ) = (z x ) = 0. Effective 
action of the gauge field and the field a x is obtained by 
integrating out the gapless spinon field z x (a = 1,2). The 
resultant effective action becomes nonlocal and therefore 
it is possible for the gauge dynamics to belong to different 
universality class from that of the usual gauge theory in 
2 + 1 dimensions. By the continuum field-theory calcula- 
tion, the effective action of the gauge field A^ is obtained 
as 
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where = d^A v — d v A^. Similarly the effective action 
of the field er T is obtained as 
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C a oc J d 3 p a(-p)-^-a(p), 
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where cr(p) is the Fourier transformed field of a x - Equa- 
tion (9) shows that fluctuations of the field a x are 
strongly suppressed at large distances. In the CP 1 model 
on the 3-dimcnsional space-time cubic lattice, a similar 
expression of the action Ca and C a is obtained by the 
hopping expansion of zf. (a = 1, 2). In the effective action 



Sa of the compact gauge field 
nonlocal terms appear 
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where the summation over closed loops T includes loops 
of an arbitrary large size, |T| is the length of T and the 
parameter 7 is estimated as 7 ~ ^ for massless zf. with 
the dimension of the lattice d = 3. We shall focus our in- 
terest on the gauge dynamics of the above nonlocal action 
which is one of the most important problems in the theo- 
retical studies on the strongly-correlated electron systems 
in these days and is still controversial. At present it is 
known that there exists only the confinement phase in 
the (2 + l)-dimensional compact U(l) gauge theory with- 
out matter couplings. 9 However phase structure is not 
clear when the U(l) gauge field couples to matter fields, 
particularly gapless matter fields 10-14 . In particular in 
Ref. 10 ' 11 , it is argued that a deconfinement phase is real- 
ized by the gapless fermion couplings. Simple mean-field 
type argument is not applicable for the nonlocal gauge 
systems and numerical studies are required. 

We shall study lattice gauge model with a nonlo- 
cal action which is related with (10) but slightly more 
tractable. The summation over T in S A (10) becomes 
(logarithmically) divergent for the configulation XI ^ = l 15 
since the massless relativistic bosons z% (a — 1, 2), which 
appear at the critical point, give divergent hopping ex- 
pansion for Up = 1. From the above discussion we shall 
consider the following (2 + 1) dimensional lattice gauge 
model with the long-range interaction in the f direction, 

Nr ~ x — 1 

Sg=9iJ2 E - U^x)W x+ll {n)Ul(x + nr)W^(n) 

(11) 
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where W x {n) — Uo{x)Uq{x + f) • • • Uq(x + nf), N T is the 
system size in the f direction and gi (i = 1,2) are cou- 
pling constants for the time and spatial directions, re- 
spectively. From action (11), it is obvious that the gauge 
model under study has nonlocal coupling in the f direc- 
tion whereas it has the usual local Maxwell-type correla- 
tion in the spatial directions. 16 Reason why we take the 
action (11) is that the f direction terms logarithmically 
divergent for £/„ = 1 and also Monte-Carlo simulations 
are easier for the model than those with full-nonlocal in- 
teraction terms. We think that studies of the model (11) 



give important insight for the full-nonlocal gauge system 
(10). More comments on this point will be given after 
showing the results of the Monte-Carlo simulations of 
the model (11). 
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FIG. 2. Fluctuation of energy as a function of gi is plotted 
for the fixed <?2 = 1.0 and 3.0. Lattice size is 12 2 x 16 and 
9 2 x 12. The results show the existence of a phase transition. 
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FIG. 3. The correlations of Polyakov lines as functions of 
their distance are plotted. The results for g\ — 0.12 and 0.14 
are shown in (a), and one for gi = 0.16 is shown in (b). Value 
of gi is fixed as <?2 = 1.0. The results show the existence of 
the confinement-deconfinement phase transition. 



We studied phase structure of the model (11) by the 
Monte-Carlo simulations of the standard Metropolis algo- 
rithm. We calculated the fluctuation of the energy E, i.e., 
((E— (E)) 2 ), as a function of the coupling g\ for fixed <?2- 
The results are shown in Fig. 2 for g2 = 1.0 and 52 = 3.0, 
and they indicate the existence of a phase transition at 
critical coupling g\c ~ 0.14 and g\c ~ 0.13, respectively. 
We varied the system size and verified that the peak gets 
sharper as the system size increases. This means that 
the transition is of second order. In order to investigate 
gauge dynamics in each phase we calculated the corre- 
lations of the Polyakov lines G{£) = (P(x)P^(x + tl)) 
where P(x) = Uo(x) ■ ■ ■ Uq(x + N t t ). 17 In Fig. 3, we show 
the results of G(£) for various values of gi and it is obvi- 
ous that the behavior of G(t) changes at g\ ~ g\c from 
]nG(£) oc £ to G(£) ~ constant for large £. Then the 
phase transition from the confinement to the deconfine- 
ment phases occurs as g\ increases. 

Numerical calculations show that the spatial coupling 
of the gauge field and larger system size enhance the de- 
confinement phase. From the above results we expect 
that the full-nonlocal gauge model (10), which results 
from the integration of the massless relativistic boson z% 
at the criticality, is in the deconfinement phase. In fact 
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the coupling constant g\ ~ 1 corresponds to (10) in which 
the damping factor balances the entropy factor of the 
paths r. 

In the deconfinement phase of (11), topologically non- 
trivial configurations are suppressed and the field-theory 
result (8) gives a qualitatively correct picture. Charges 
interacting through have the potential V(r) oc 1/r 
where r is the spatial distance between the two charges. 

Let us comment on the effects of the Berry phase. Since 
Cb is neglected in (11), one may doubt the deconfinement 
phase transition observed above. However, as the Berry 
phase generates extra phases for topologically nontrivial 
configurations in the path integral, the Berry phase en- 
hances the deconfinement. In fact without these extra 
phases, all instanton configurations contribute additively 
to disorder the gauge system. Thus the existence of the 
deconfinement phase in the gauge system (11) guaran- 
tees its existence even in the presence of the Berry phase. 
Similar argument was used for the deconfinement transi- 
tion at finite T 18 , which is established at present. 

We summarize the phase structure of the original spin 
model. In the region a < etc (/cflf < fc), the low-energy 
excitations are the massless spin waves whereas in the 
region a > ac {f e S > fc), they are s = 1 excitations 
n x . On the critical point a — ac (,f c g = fc), the gauge 
dynamics is in the "Coulomb" phase and the low-energy 
excitations are the s = 1/2 bosonic spinons z x (a = 1,2) 
which are interacting with each other by the potential 
1/r. The spin correlation function decays algebraically 
both in the Neel state and at the criticality but exponent 
is different. In the Neel state, the spin operator is given 
as S x = \z x az x <~ (z x )az x ~ z x whereas at the criticality 
n x , the bilinear of z% and z x ^ (a = 1,2). Phase structure 
of the nonuniform AF Heisenberg model is schematically 
shown in Fig. 4. 



theories will be reported in near future. 20 
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FIG. 4. Schematic phase diagram of the Heisenberg model. 
Phases A,B and C are the Higgs (Neel), Coulomb (critical) 
and confinement (dimer) phases, respectively. 



Recently a similar phase transition from the Neel to 
the dimer states was discussed 19 . There they conclude 
that instanton effects are irrelevant at the critical point. 
Our numerical investigation is consistent with their result 
but our study shows the long-range interactions of the 
gauge field, which appear as a result of the coupling to 
the massless boson z x , play an essentially important role. 
Results of more detailed studies on the long-range gauge 
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